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Primordial non-Gaussianities produced by features in the potential of 
single slow-roll inflationary models 

by Alexander Gallego 


We study the effects of a general type of features of the inflaton potential on the 
spectrum and bispectrum of primordial curvature perturbations. These features 
correspond to a discontinuity in the n-th order derivative of the potential which are 
dumped exponentially away from the value of the held where the feature happens. 
Interestingly we And that different values of the amplitude and of the order of 
the feature can lead to the same effects on the power spectrum on both large and 
short scales, and on the bispectrum at small scales. Only taking into account the 
bispectrum at large scales it is possible to resolve this degeneracy. We provide 
fully numerical calculations and analytical approximations for the spectrum and 
the bispectrum, which are in good agreement with each other. The analytical 
approximation allows to to determine the class of features which give the same 
spectrum and can only be distinguished with the bispectrum. 
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Chapter 1 


Introduction 


In the last few decades the outstanding advances in observational cosmology have 
set up the Standard Model of Cosmology [1-3]. There is a signihcant set of observa¬ 
tional data from experiments like the Cosmic Background Explorer (COBE) Satel¬ 
lite, the Sloan Digital Sky Survey (SDSS), the Wilkinson Microwave Anisotropy 
Probe (WMAP), and the Planck mission [2, 3]. There are also many other ground- 
based and sub-orbital experiments [4, 5]. Recent experiments reveal a universe that 
is 13.81Gyrs old, and made of 4.9% Baryonic Matter. The rest is made of Dark 
Matter, 26.8%, and 68.3% of a unknown type of energy, called Dark Enery [1-3, 6]. 

According to the Cosmological principle, the universe is isotropic and homogeneous [7, 
8]. However, results from the CfA survey, and the SDSS, reveled an inhomoge¬ 
neous distribution of the galaxies in the observable universe. This showed that 
the universe is only homogeneous at large scales (at about lOOMpc). Moreover, 
results from COBE, WMAP, and Planck showed that our universe is not perfectly 
isotropic [3, 9, 10]. 

The Cosmic Microwave Background (CMB) radiation is the radiation coming from 
the last scattering surface emitted when the universe was only about 378, 000 
years old. At this time, protons and electrons combined to form neutral light 
atoms (recombination) and then photons started to travel freely through space 
(decoupling) [2, 10]^. Although this radiation is extremely isotropic, it is still 
possible to observe small fluctuations in the temperature map, determined by 
AT /T ~ 10“® [2, 9-11]. The universe was very young at the recombination epoch. 

^Recombination and decoupling occurred at different times [7, 8]. 
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Thus, a detailed measurement and subsequent study of this radiation can gives us 
a very valuable information about the physics of the early universe [11]. 

In a perfectly homogeneous and isotropic universe, the process of formation of 
stars, and development of more complex structures (due to gravitational instabil¬ 
ity) would never have begun [7, 12, 13]. Thus, any model trying to explain the 
evolution of the universe must include some level of inhomogeneity to account for 
the formation of Large-Scale Structures and anisotropies in the CMB. For instance 
WMAP and Planck data are very well htted by the ACDM model [2, 3]. And many 
different inflationary models can account for the formation of anisotropies in the 
universe [11, 14]. 

Cosmological inflation is a period of the universe’s evolution during which the scale 
factor is accelerating [12, 15, 16]. But instead of being a replacement for the Big 
Bang, inflation is actually more like an odd-on: after the primordial inflationary 
period, the universe returns to a non-accelerating phase in which all the success 
of the Big Bang model is recovered [12]. 

During inflation, space-time itself fluctuates quantum mechanically about a back¬ 
ground in an accelerated expanding phase. These microscopic fluctuations were 
spread out to macroscopic scales, where they eventually become small classical 
fluctuations in the space-time. From then on, slightly over-dense regions started 
the process of gravitational collapse, forming the stars and galaxies [6-8, 12]. 

For a long time, observations of the fluctuations in the CMB and the Large-Scale 
Structures have been focused mainly on the Gaussian contribution as measured by 
the two-point correlations of the density fluctuations (or by its Fourier transform, 
known as the power spectrum) [1, 12, 17]. It has been precisely the study of the 
power spectrum which has settled the numerical value of many of the cosmological 
parameters in the ACDM model [1]. 

The main justihcation to go beyond the study of the Gaussian primordial fluc¬ 
tuations is discriminate between models of inflation which can generate the same 
power spectrum [11, 18, 19]. Thus, in order to distinguish these models we study 
the generation of fluctuations in slow-roll inflationary models, and then com¬ 
pare these results with the non-Gaussianity on the correlation functions obtained 
experimentally[l-3, 9-11, 13]. 



Chapter 2 


The Standard ACDM Model of 
Cosmology 


When we look at the sky we see a nearly isotropic universe. This means that 
the universe looks almost the same in all directions. Observations also indicate 
that at large-scales the universe is homogeneous [20, 21], which means that if we 
move from one random point to another in the universe then everything looks 
the same. The cosmological principle is based on these empirical observations. It 
states that at large-scales the visible parts of our universe are homogeneous and 
isotropic [22, 23]. 

In 1916, Einstein published his theory of General Relativity (GR), which related 
the geometry of space-time with the energy-matter content of the universe. A 
few years later, Alexander Friedmann (1922) and Georges Lemaitre (1927) in¬ 
dependently applied the cosmological principle to GR. They solved the Einstein 
equations by assuming a homogeneous and isotropic spacetime. The resulting 
equations revealed the possibility of an expanding universe. In a 1927 article 
Lemaitre even suggested an estimated value of the rate of expansion. It was only 
two years later in 1929 that Edwin Hubble conhrmed the existence of this expan¬ 
sion and determined its rate, now called the Hubble rate parameter H. It was the 
birth of modern cosmology: for the hrst time in history we had a compelling and 
testable theory of the universe. 

At the beginning of this chapter we assume that the universe is described by the 
Big Bang model of cosmology [24, 25]. In Section 2.2 we review a series of classical 
observational tests which conhrm the predictions of the model. Then in Section 
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2.3 we point out some shortcomings of the Big Bang that will lead us beyond the 
standard Big Bang model. The cosmological inflationary theory is introduced in 
Section 2.3 in order to solve the problems that the Big Bang model is insufficient 
for. The current standard model of cosmology called the ACDM model [2, 3] is 
described in Section 2.5. 


2.1 The Big Bang Model 


The Big Bang model which rests on two theoretical principles: the cosmological 
principle and General Relativity as the correct theory of gravity on cosmological 
scales. In this model the energy-matter content of the universe is divided into 
four primary components [22, 23]: i) radiation, which is composed of relativistic 
or nearly relativistic particles such as photons and neutrinos, ii) “ordinary mat¬ 
ter” composed of protons, neutrons, and electrons. It is generally called baryonic 
matter, hi) cold dark matter (CDM) which refers to exotic non-baryonic matter 
that interacts only weakly with ordinary matter, and iv) dark energy (DE), char¬ 
acterized by a large, negative pressure. This is the only form of matter that can 
accelerate the expansion of the universe. The hrst three components are commonly 
united in just one form simply called the matter component [3]. 

The Einstein held equations describing the dynamics of the universe with metric 
tensor g^i, are given by 


( 2 . 1 ) 






where we have taken the speed of light c = 1. and are the Riemann and 
energy-momentum tensors, respectively. TZ is the Ricci scalar, and A and G are 
the cosmological and gravitational constants, respectively. The metric signature 
that we adopt is (—1,1,1,!). When the cosmological principle is assumed the 


spacetime takes the form of the Friedmann-Lemaitre-Robertson-Walker (FLRW) 
spacetime, with metric 


= g^i.dx^dx'' = —dC 



where a{t) is the scale factor [11, 26]. The constant of curvature is denoted by 
K with three possible set of values K < 0, J’C = 0, and K > 0, corresponding 
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to open, flat, and closed universes [26], respectively. The Einstein held equations 
give 


Friedmann equation 

/aV StiG K 
\ a / 3 ^ a^ ’ 

(2.3) 

Acceleration equation 

- =-^(p + 3P), 

a 3 

(2.4) 


where dots indicate derivatives with respect to time. The Hubble rate parameter 
which measures how rapidly the scale factor changes is defined as if = a/a. From 
the energy conservation condition = 0 we can derive a third equation which 

will be useful later 

Fluid equation ^ + 3- (p + P) = 0, (2.5) 

dt a 

where p and P are the energy density and pressure in the universe, respectively. 
Here p is the sum of several different components of species: radiation, baryonic 
matter, CDM, and DE. The fluid equation expresses the energy conservation for 
the fluid as the universe expands. 

From the equations of motion it is clear that the rate of expansion H is related to 
the energy-matter content of the universe p and to the geometry of spacetime K. 

The fluid equation Eq. (2.5) can be written as 

2(po=>) + p2o=> = 0. (2.6) 

The comoving volume H = a^ is the volume of a region expanding together with 
the cosmic fluid, thus the energy in a comoving volume is given by P = pa^. Then 
Eq. (2.5) takes the simple form 


dU + PdV = 0. (2.7) 

For a fluid in equilibrium the first law of thermodynamics states 

dU + PdV = TdS, (2.8) 

where T is the temperature and S the entropy [8]. An adiabatic process is one 
in which dS = 0. This is precisely the case of Eq. (2.7) which shows that a 
homogeneous and isotropic universe with a perfect fluid expands adiabatically 





Chapter 2. The Standard ACDM Model of Cosmology 


6 


[27]. In case that the perfect fluid obeys the barotropic relation 

P = wp, 


( 2 . 9 ) 


Eq. (2.6) is expressed as 




( 2 . 10 ) 


In this chapter we adopt the usual convention of writing the values at the present 
time with a subscript 0. The solution is given in terms of the present value of the 
density po 


P 


Poa 


3(ii)+l) 


( 2 . 11 ) 


For radiation we have w = 1/3, for matter, which have effectively zero pressure, 
ta = 0, and for a cosmological constant w = —1. Thus 


P-y = Po,'ya^, ( 2 . 12 ) 

Pm ~ P0,mOi , (2.13) 


where pv is constant. As the scale factor increases, the density of radiation will 
decrease faster than for matter. The density of a cosmological constant does not 
change as the universe expands. Recent observations [2, 3] show that today matter 
and a cosmological constant dominate the energy content in the universe, but since 
p^ oc as a decreases as time goes back, p increases, thus radiation was the 
dominant content in the early universe. There is a special value of p for which 
A = 0, it is known as the critical density 

. , 3A2(t) 

pM . ( 2 . 14 ) 


An useful dimensionless quantity called the density parameter is dehned as 


fi(f) 


pif) 

Pc{t) ■ 


(2.15) 


With this definition the Friedmann equation Eq. (2.3) can be written as 


^tot I5 


( 2 . 16 ) 


where ^Itot = ^m + is the sum of all the energy-matter components and 12^ = 
—K/{aH)‘^ is the density parameter associated with the curvature term. Planck 
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reported the values for him and parameters by htting the model to the data 
[3, 25]. The best ht parameter values for the density parameters are Qm = 0.314± 
0.020 and Ha = 0.686 ± 0.020. From this results and from Eq. (2.16) it is safe to 
conclude that the universe is nearly flat [3]. Thus, from now on, we consider a flat 
universe K = 0 with spacetime metric 


= —dt^ + a{tYdx‘^, 


( 2 . 17 ) 


where x are the comoving spatial coordinates and the fundamental distance mea¬ 
sure is the distance on this comoving grid. For example, an important comoving 
distance is the distance light could have traveled since a time tj. In a time dt, light 
travels a physical distance dt = adx so the total comoving distance light could 
have traveled, which is called comoving horizon is 



( 2 . 18 ) 


It is called the comoving horizon because no information could have propagate 
further than r in the comoving grid since tj. If we now consider = 0 as the 
beginning of time then regions in the universe separated by a distance greater 
than r are not causally connected. In the next chapters we will use r, also called 
the conformal time, to study the evolution of perturbations. 

In the case of a flat universe the evolution of the scale factor in terms of time 
is easy to hnd. Substituting the density (2.11) into Eq. (2.3), the Friedmann 
equation can be written 



( 2 . 19 ) 


which can be further written as 



( 2 . 20 ) 
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The term in the square root is constant in time. If we choose the scale factor at 
the present time fo as = 1, the solutions for the different components are 


^ id 

/ ( X 2/3 

^ u 


a Ait) = exp 


SvrpA 


it - to) 


( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 


As can be seen from the last equation a constant energy density gives a very rapid 
expansion. And the Hubble constant for each component is 


H,it) 

Hmit) 

HAit) 


1 

2 



= constant. 


(2.24) 

(2.25) 

(2.26) 


If we suppose that the universe have been matter-dominated by most of its time, 
then by measuring the Hubble rate today Hq = (2/3)fQ^ and the age of the 
universe to separately we could have a powerful test for the Big Bang model. In 
the following section we will confront the theory of an expanding universe with the 
astronomical observations. After this tests we will see the need of extending the 
Big Bang theory by adding cosmological inflation, the new model will be called 
the ACDM model [3]. 


2.2 Observational Tests 

Nowadays we have good evidence that the universe is indeed expanding. This 
means that early in its history the distance between galaxies were smaller than 
they are today. It was then an extremely hot and dense state which began to 
expand rapidly. This early development of the universe was known as the Hot Big- 
Bang. The success of the Big Bang theory rests on three observational pillars: i) 
expansion of the universe, ii) nucleosynthesis of light elements, and hi) origin of the 
cosmic microwave background (CMB) radiation. We review now this measurable 
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signatures which strongly support the notion that the universe evolved from a hot 
dense, homogeneous, and isotropic gas as proposed by the Big Bang model. 


Expansion of the universe According to the Big Bang the universe is ex¬ 
panding in all directions, as it can be seen in the equations of motion from GR 
[Eqs. (2.3)-(2.4)]. The Hubble diagram Fig. 2.1 is the most direct evidence of 
this prediction [3, 28]. The Hubble’s law states that the recessional velocity v of 
a galaxy is proportional to its distance d from us, v = Hod. Current values of the 
Hubble constant are parameterized by the little h as [3, 7] 

Hq = lOOh km s“^Mpc“^, (2.27) 


where Mpc denotes the astronomical length scale megaparsec equal to 3.086 x 10^^ 
m. From Planck results [3] the present value of h is h = 0.673 ± 0.012. As the 
galaxies are receding from us their emitted wavelength \emit should be stretched 
out so that the observed wavelength Xobs is larger than Xemit- A measure of this 
stretching factor is called the redshift 2 ; dehned by 


z = 


'06s 


A 


emit 


^emit 


ao — a 
a 


(2.28) 


and since oq = 1, the farthest objects are at larger redshift. Moreover, current 
experiments show that the expansion rate of the universe is accelerating [29, 30]. 
The discovery of cosmic acceleration has indicated the presence of DE in the 
universe [2, 3], since, as we saw above in Eq. (2.23), a cosmological constant can 
accelerate the expansion rate of the universe. It also has stimulated new physics 
beyond the standard model from modihcations to GR as possible explanations of 
its presence [19]. 


Big-bang nucleosynthesis The Big-bang nucleosynthesis (BBN) is a nonequi- 
librium process that took place just a few minutes after the Big-Bang when nu¬ 
cleons were synthesised into the light elements such as D, ^H, ^He, and ^He. To 
study the BBN we need some basic assumptions from GR and the standard model 
of particle physics. Less than one second after the Big-Bang the universe was very 
hot with a temperature larger than lO^'^K. It was a rapidly expanding plasma, 
with most of its energy in the form of radiation and high energy relativistic par¬ 
ticles, such as a photons and neutrinos. This high radiation made impossible the 
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Distance (Mpc) 

Figure 2.1: Hubble diagram from the Hubble Space Telescope Key Project. 
Top: Hubble diagram of distance vs. velocity for secondary distance indicators 
calibrated by Cepheids using five different measures of distance. Bottom: Value 
of Hq as a function of distance. Image credit: Freedman et al. [28] 


formation of any atom or nucleus for if it were produced, it would be soon de¬ 
stroyed by these relativistic particles. As the universe expanded the temperature 
dropped well below the binding energies of typical nuclei MeV/fcg giving rise to 
the light elements [25, 31]. These predictions of the BBN are in good agreement 
with current observations [25]. For instance in the case of primordial abundance 
of deuterium these observations correspond to probes of clouds at high redshift on 
the line of sight of distant quasars. While ^He is determined from observations in 
ionized hydrogen regions of compact blue galaxies [25]. 

In the Big Bang model one of the most important parameters is the baryon to 
photon ratio, r] = n^/n^ [25, 31] where the photon number density is determined 
from the CMB temperature and is related to the baryonic density fife [25]. Fig. 
2.2 shows the abundances of ^He, D, ^He and ^Li from the stellar and extragalactic 
observational experiments (as the ones mentioned above) and from the WMAP and 
Planck missions [25]. From these and other observations the successful predictions 
of BBN make it a cornerstone of the Big Bang cosmology [3, 25, 31]. 
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Figure 2.2: Abundances of '^He, D, ^He and ^Li (blue lines) as a function of 
(top) or g (bottom) [25]. The observational abundances are represented by 
the green horizontal areas. The WMAP and Planck baryonic densities results 
are represented by the dotted-black and solid-yellow lines, respectively. Image 
credit; Coc et al. [25]. 


Cosmic microwave background Radiation In addition to the assumption 
of a homogeneous and isotropic universe, we also consider a universe whose rate 
of expansion H is less than the interaction rate per particle V = nav [8], where 
n is the number density of the target interacting particles, a is the cross section, 
and V is the average relative velocity. This last assumption tells us that after the 
Big-Bang particles were in thermal equilibrium [8, 32]. Which is not surprising 
since homogeneity and isotropy imply no temperature gradients and hence no heat 
flow. Photons then remain tightly coupled to electrons via Compton scattering. 
Any hydrogen atom produced was immediately ionized due to the low baryon to 
photon ratio rj. As the temperature dropped below leV, when the universe was 
roughly 378,000 years old, protons and electrons combined to form neutral light 
atoms. This period is known as the recombination epoch. Decoupling occurred 
shortly afterwards, when photons started to travel freely through space since their 
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Figure 2.3: Brightness of the cosmic microwave background radiation as a 
function of wavelength (top) or frequency (bottom) from Differential Microwave 
Radiometers (DMR), Far Infrared Absolute Spectrophotometer (FIRAS) and 
other ground-based and sub-orbital experiments. The theoretical curve for a 
blackbody spectrum is represented here by the dashed line. Image credit: G. 
Smoot [34]. 


energy was not enough to continue been scattered by electrons. These photons 
constitute the relic radiation coming from the last scattering surface that we ob¬ 
serve today. And, since before decoupling photons were at thermal equilibrium, 
this radiation should have a blackbody spectrum today. Thus when we observe 
them, we see a uniform radiation coming from all directions in the sky. The 
best-stimate CMB photon temperature today is Tq = 2.72548 ± 0.00057K [3, 33]. 
Fig. 2.3 shows the spectrum of the CMB radiation for various experiments (See 
Ref. [34] and references therein). The dashed line corresponds to the theoretical 
curve of a blackbody spectrum. The agreement between theory and experiments 
is remarkable and is one of the most powerful predictions of the Big Bang. 


2.3 Problems with the Big Bang Model 


The Big Bang theory has encountered remarkable successes, however it is not com¬ 
plete, and a few unsatisfactory facts remain unexplained. It seems that the model 
necessitates a larger theoretical framework. We review now the main problems: i) 
flatness problem, ii) horizon problem, and hi) primordial perturbations. 
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The flatness problem In the Big Bang model, the scale factor behaves like 
a (X C with q < 1, see Eqs. (2.21) and (2.22). Rewriting Eq. (2.16) in the form 


l^tot 1 | ^ 2^2 , 


(2.29) 


it can be seen that \C,tot — 1| is an increasing function of time in either case. That 
means that the flat geometry is a unstable situation for the universe: any small 
departure from this stringent initial condition would quickly turn the universe into 
open or closed. And considering the current value of Qtot ~ 1 from Planck [2, 3] 
it requires an extreme flne-tuning in the early universe, something that the Big 
Bang paradigm cannot explain. 

The horizon problem To understand the horizon problem consider the max¬ 
imal physical distance that can be covered by a light ray dhor = or since the 
beginning of time ti = 0. From Eq. (2.18) this distance is given by (remember 
that a oc C) 



(2.30) 


which is called the (particle) horizon [35] and it is of the order of the so-called 
Hubble radius We thus have that for a flnite time dhor is flnite, meaning that 
the scales at which particles or perturbations were connected at early times had to 
be small. How could then different regions in the universe, which are separated by 
large distances, have nearly the same temperature and other physical properties 
if they could not have been causally connected? 

Primordial perturbations The Big Bang assumes a uniform, isotropic uni¬ 
verse from its beginning. But this uniformity would not allow the formation of 
galaxies, or clusters of galaxies in the universe. We know that gravitational insta¬ 
bility is the responsible for the structure in the universe. After the Big Bang, as 
the temperature continued to descend due to the expansion, the energy density of 
the universe began to be dominated by matter instead of radiation. The primor¬ 
dial perturbations in the matter density accumulated matter from the surrounding 
areas. In this way initially overdense regions began to grow, independently of how 
small this overdensity was. Correspondingly, the slightly underdense regions grew 
more underdense. Hence we need an explanation of how this primordial fluctu¬ 
ations were created in order to give rise to underdense/overdense regions in the 
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Figure 2.4: Comparison of the three satellite’s measures of the CMB temper¬ 
ature. The image shows tiny small variations (anisotropies) in the temperature 
readings, specially from WMAP and Planck missions. Hotter (colder) regions of 
the sky are represented by red (blue). Image credit: NASA/JPL-Caltech/ESA. 

universe. A detailed search for these primordial perturbations started with the 
COBE satellite, then with the WMAP, and culminated in the recent Planck mis¬ 
sion as can be seen in Fig. 2.4. Hotter (colder) regions of the sky are represented by 
red (blue). The temperature fluctuations correspond to regions of slightly different 
densities at early times, representing the seeds of all future structure: clusters of 
galaxies, galaxies and stars [9, 10]. In the COBE mission the resolution to measure 
the fluctuations were of about 7°. While WMAP managed to go down to about 
0.5°. The amazing results from Planck come from a resolution of 0.1° which allows 
it to measure temperature fluctuations down to a 10“®K [36]. The Planck map of 
the entire sky is shown in Fig. 2.5. 

In the next section we will see how the inflationary cosmology theory, when added 
to the Big Bang, can solved the shortcomings mentioned above. 


2.4 Solutions to the Big Bang problems from In¬ 
flation 

In order to solve the problems of the Big Bang it is necessary to introduce cos¬ 
mological inflation [3, 10, 15]. In the next chapter we will propose a precise and 
quantitative theory of inflation characterized by a scalar held as the energy-matter 
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Figure 2.5: All-sky image of the Cosmic microwave background as seen by 
Planck. Hotter (colder) regions of the sky are represented by red (blue). Image 
credit: ESA and the Planck Collaboration. 


component of the universe before the Big Bang. For the moment we will only con¬ 
sider its qualitative dehnition and implications. 

Inflation is any period of the universe’s evolution during which the scale factor is 
accelerating 

INFLATION h > 0. (2.31) 

We will see that from this simple dehnition the unsolved problems of the Big Bang 
theory can be relieved. 


The flatness problem This is the simplest problem that can be solved. From 
the inflationary condition Eq. (2.31) we have that 

a > 0 ^d > 0 (2.32) 

LLL LLL 

which indicates that during inflation aH is an increasing function of time. Thus 
from Eq. (2.29) we conclude that as aH increases \ fltot — 1| has to decrease, i.e., 
Qtot ~ 1 in concordance with Planck results [3]. 


The horizon problem The Hubble radius is the distance over which particles 
can travel in the course of one expansion time [7]. Another useful quantity is the 
comoving Hubble radius H~^/a = {aH)~^. If particles are separated by distances 
greater than {aH)~^ they cannot communicate with each other today, but it is 
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possible that they were in causal contact early on [35]. As opposed to the comov- 
ing horizon r in which particles that are separated by a distance larger than r 
could have never communicated, since it is the maximal comoving distance travel 
by light. A solution to the horizon problem comes from a shrinkage of the comov¬ 
ing Hubble radius so that particles that were initially in causal contact with one 
another early on can now no longer communicate. We thus would need {aH)~^ 
to decrease, then aH must increase. But this is simply the inflationary condition 
met above since 

^iaH) > 0 -^a > 0 a > 0. (2.33) 

dV ^ dt ^ ^ 

Primordial perturbations The quantum fluctuations of the inflaton held gen¬ 
erates the primordial perturbations responsible for the formation of structures. 
According to this, the universe should exhibit anisotropies from its early times. 
Fig. 2.6 shows the anisotropy spectrum (power spectrum) of the CMB temper¬ 
ature from Planck [36]. Planck has now measured the hrst seven acoustic peaks 
to high precision that are well ht by inhation and the standard Big Bang model 
which we will called ACDM [2, 3]. 


2.5 The Standard ACDM Model of Cosmology 

The remarkable agreement between experiments and the Big Bang plus the in¬ 
hationary theory have set up the standard model of cosmology. From now on 
we will simply call it the ACDM model [2, 3]. The model is based upon an 
expanding universe whose dynamics are governed by General Relativity with a 
spatially-hat metric and whose constituents are dominated by CDM and, recently, 
by a cosmological constant A. The primordial seeds of structure formation gener¬ 
ated during inhation are adiabatic and have a gaussian distribution with a nearly 
scale-invariant spectrum [3]. This model is described by only six parameters. 
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Figure 2.6: The temperature angular power spectrum of the primary CMB 
from Planck [36]. The vertical scale is P; = l{l + 1)Ci/2tt vs the angular scale 
(top) or the multipole-moment (bottom). Cosmic variance is represented by 
the shaded area and is also included in the error bars on individual points. The 
horizontal axis is logarithmic up to ^ = 50, and linear beyond. The green line 
is the theoretical prediction from the ACDM model. Image credit: Ade et al. 
Planck 2013 results [36]. 










Chapter 3 


Cosmological Perturbations 


The cosmological inflationary scenario provides a possible mechanism for gener¬ 
ating the spatial and matter variations in the very early universe. These are the 
seed fluctuations for the observed structures and CMB anisotropies, hence it is 
the finest and leading theory explaining the early universe [1, 12, 14-16]. 


3.1 Inflation 

From the acceleration equation, Eq. (2.4), we can see that in order to fulfill the 
inflation condition a > 0 we must need that 

P < (3.1) 

and since the energy density is always positive, the pressure must be negative. A 
negative pressure can be modeled by a cosmological constant A with tc = — 1 as we 
saw before, Eq. (2.9). Also from Eq. (2.23) we see that A induces an accelerated 
expansion. However, inflation must come to an end after some amount of time [10], 
thus A would have to decay into the “conventional” particles [12]. And postulating 
that a cosmological constant can decay away, converting its energy into another 
forms is not a viable way to bring inflation to an end [12]. Besides, we do not want 
to spoil the success of the Big Bang model. So in typical models, inflation occurs 
when the universe is extremely young, before nucleosynthesis, perhaps around 
10“^^s. Hence, at this large energy scale of around lO^^GeV we must consider the 
fundamental interactions of particles. 
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3.1.1 The inflaton 


As the universe expands it cools down. It is known that when a physical system is 
heated or cooled the properties of the system can change dramatically, leading to 
phase transitions [12]. In fact, an example of a phase transition could have been 
when quarks hrst condensed to form hadrons [8]. In fundamental particle physics 
a phase transition can be controlled by a form of matter known as a scalar held. 
Thus, in the simplest inhationary universe, a scalar held would decay away after 
the phase transition, ending inhation and bringing a large amount of expansion in 
order to solved the hatness and horizon problems. And since we are considering a 
held at very high energies, we know from quantum mechanics that the held have 
quantum huctuations in spacetime. Moreover, spacetime itself huctuates quantum 
mechanically about a background that is expanding at an accelerating rate. These 
huctuations are then responsible for the generation of matter and spatial variations 
in the very early universe [12, 13]. 


We consider the simplest inhationary model, composed of a single scalar held cj) 
called the inflaton^. We hrst write down the energy-momentum tensor for 0 and 
then investigate if this kind of held could have a negative pressnre P and then 
if it is actnally responsible for an accelerated expansion. The energy-momentum 
tensor can be written [12] as 


Tfxu — d^(j)dy(j) g^j^p 




(3.2) 


where V{(j)) is the potential energy for the held. The energy density p is the 
time-time component 

p = \i>'‘+yw- ( 3 . 3 ) 

The hrst term is the kinetic energy density. The pressure P is T\ (no sum over 
i), thus we have 

p = iy - V(,l,). (3.4) 

Thus, for a scalar held to have negative pressure it is necessary that its kinetic 
energy be less than the potential. The analogy of a scalar held trapped in a 
false vacuum to achieved a negative pressure was used in the initial formulation 
of inhation by A. Guth [15]. Bnt it was soon realized that such scenario was 
unfeasible [37]. To avoid this and several other problems as a conseqnence of a 


Mhe case of multi-field inflation can be studied in Refs. [11, 35]. 
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trapped field the infiaton was considered slowly rolling toward its true vacuum [7]. 
The small change in its kinetic energy implies 0 is almost constant and then its 
total energy is all potential V, which remains constant with time. A spacetime 
with a constant positive vacuum energy density Vq is called a de Sitter space. We 
have already seen the effect of constant energy density on the evolution of the 
universe, Eq. (2.23). In the case of a universe dominated by a scalar field density 
p we have that 

a ISttGp 

- = H = \ -cs constant, (3.5) 

a V 3 

hence 

a{t) ^ a{ti)e^^^~^^\ (3.6) 

where is a beginning time for inflation and a(fj) the initial value of the scale 
factor. Thus the infiaton cj) can be responsible for an accelerated expansion. 

We first study the zero-mode background evolution of spacetime and the infiaton. 
Then we move onto perturbations in spacetime and the scalar field (j). 


3.1.2 The background evolution 


We could in principle determine the dynamics of 0 in the background from the Ein¬ 
stein equations by simply replacing the energy density and pressure of the infiaton 
[Eqs. (3.3) and (3.4)] into Eqs. (2.3)-(2.4). In this chapter we follow a different 
approach though, since we should also consider later the quantum dynamics of the 
infiaton and the spacetime itself under perturbations. And this approach is more 
conveniently achieved by the ADM formalism in which the starting point for the 
formulation is the Lagrangian [11, 38]. 


The dynamics of the spacetime and the scalar field 0 driving inflation are governed 
by the action [12, 35] 


_ Mpi 
~ ^ 


d^XyJ—gR 





(3.7) 


where Mpi = (SvrG)”^/^ is the reduced Planck mass. The first term is the Hilbert- 
Einstein action. The second term represents the minimal coupling of the infiaton 
to gravity through the background metric tensor [11, 12, 35] 


ds^ = g^i/dx^dx'" 


dt^ -I- a{tYdx‘^, 


(3.8) 
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where x are the comoving spatial coordinates. The classical equations of motion 
are given by 

/ • \ 2 1 /I \ 

/ n \ \ / I . - \ 

(3.9) 


3MJ, 


ly +!/(.).) 


(j) + SHcj) + d^V — 0, 
^2 

H = - 


2M|,’ 


(3.10) 

(3.11) 


where dots and indicate derivatives with respect to time and the scalar held, 
respectively. The last equation is redundant since it can be derived from the hrst 
two; it is called the the continuity equation. 


The amount of inhation is quantihed by the ratio of the scale factor at the hnal 
time to its value at some initial time. Normally, this ratio is a large number, thus 
it is customary to take the logarithm to give the number of e-foldings 


N{t) = In 


a{ti) 


(3.12) 


It is also called the number of Hubble times, since N = H{t')dt'. In order to 
obtain at least N ~ 60 [11, 12], we need to impose that the Hubble parameter H 
does not change much within a Hubble time H~^. This requisite gives rise to the 
hrst slow-roll condition [11, 13, 35] 




(3.13) 


The second slow-roll condition is given by the requirement that e does not change 
much within a Hubble time [11] 



0 ( 1 ). 


(3.14) 


From the dehnition equations it can be seen that the slow-roll parameters are 
dimensionless. By applying the background solution Eq. (3.11) the slow-roll 
parameter e can also be written as 


e 


0^ 


« 0 ( 1 ), 


(3.15) 


where we can see that the kinetic energy of the scalar held is neglected, thus 
the energy driving inhation is dominated almost entirely by the potential V [11, 
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35]. From the discussion of section 3.1.1 this result was expected. The Hubble 
parameter can be approximated then by 




3Mj, 


H. 


(3.16) 


The second condition can also be rewritten as 


V 



« 0 ( 1 ), 


(3.17) 


SO (j) is negligible and Eq. (3.10) be expressed by the attractor equation [11] 


+ d^V ^ 0 . 


(3.18) 


The slow-roll parameter thus described impose restrictions on the rolling veloctity 
of the scalar held 0. They are also called the Hubble how function (HFF) slow-roll 
parameters and can be interpreted as measure of the deviation from a purely de 
Sitter background [10]. For instance if the scalar held is exactly constant then 
e = 0, and the spacetime is a purely de Sitter space given an exact exponential 
expansion. The hierarchy of these parameters is dehned as ei = —H/H‘^,ei+i = 
ei/{Hei), for i < 1. An equivalent dehnition of the slow-roll parameter is given by 
[ 10 ] 


ey = 


_ Mj, fd^V 


1 / 


« 0 ( 1 ), 


„ ^ « 0 ( 1 ). 


(3.19) 

(3.20) 


By using (3.18) they are related to the HFF slow-roll parameters by [10] 


e = ey, r^y =-2r7y-F 4ey. 


(3.21) 


Although the various dehnitions of the slow-roll parameters are equivalent, the 
dehnitions Eqs. (3.13) and (3.14) are more general [11]. 
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3.2 non-Gaussianity 

An important prediction of inflationary cosmology is that there should be small 
departures from the large-scale homogeneity observed in the universe. Inflation 
also predicts that these perturbations have a characteristic spectrum and bispec¬ 
trum. During inflation, space-time itself fluctuates quantum mechanically about a 
background in an accelerated expanding phase. These are the seed fluctuations for 
the observed structures and CMB anisotropies. The microscopic fluctuations were 
spread out to macroscopic scales, where they eventually become small classical 
fluctuations in space-time. From then on, slightly over-dense regions started the 
process of gravitational collapse, forming the stars and galaxies [ 6 - 8 , 12 ]. 

To understand the primordial fluctuations of a quantum held we need to dehne 
a quantity which brings into mutual relation its huctuations in different places. 
That quantity is given by the n-point correlation functions [13, 39] 

(D(f) I C(^i, t)C{x2, t).. .({xn, t) I n{t )), (3.22) 

where these functions represent the extent to which huctuations in diherent places 
are correlated with each other in the ground state of the interacting theory |D(t)). 
In case of the early universe, the huctuations are small compared with the back¬ 
ground, then the correlation functions can be evaluated by perturbation theory 
[13, 39]. The dynamics governing the held C{x,t) is determined by the action S 
with the associated Lagrange density C [C]. Expanding S' as a series in powers of 
the perturbation held we have 

S = I d‘x {£<») 1C] + £<‘) Id + Id + Id + ... (3.23) 

A theory is called Gaussian if all the higher order interaction Lagrange beyond 
the quadratic are zero [12, 13, 40], i.e. if [(C] = 0 for i > 2. We have that 
higher oder interactions beyond the cubic one can be written in terms of powers 
of the lower order correlators [13], thus in the case of a Gaussian distribution we 
have that 

S= /'a{£I“>[C1 + £'"’[C1}. (3.24) 

is an exact expression. Any departure from Gaussianity is called non-Gaussianity 
[1, 12, 13], i.e. if [C] 7 ^ 0 for i > 2. In Eq. (3.24) the term [((■] vanishes 
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since the background evolution is always a solution to the equations of motion 
[13, 17], 

In many single-held slow-roll models the inhaton is weakly coupled to gravity, 
hence the predicted non-Gaussianity is very small and it would be very hard to 
detect [1, 6, 35]. There is, however, a mechanism for producing a large, and 
detectable amount of non-Gaussianity, which corresponds to the violation of any 
of the following conditions [1]: 

1. Single field. In the generation of the primordial seed huctuations, there 
was only the action of one quantum held, (j). This held was responsible for 
driving inhation (there is no consideration to multi-helds). 

2. Slow-roll. During inhation, the evolution of the scalar held was very slow 
compared to the Hubble time 

3. Canonical kinetic energy. The kinetic energy of the quantum held is 
such that the speed of propagation of huctuations is equal to the speed of 
light. 

4. Initial vacuum state. The quantum held was in the “Bunch-Davies” 
vacuum state [11], right before the generation of the quantum huctuations 
during the inhationary phase. 

In chapter 4 we will study the introduction of a feature in the potential of the 
inhaton which leads to an enhancement of the the slow-roll parameter rj. In this 
case e is small but g is large thus, due to this violation of the slow-roll condition, 
we expect a large contribution to the three-point correlation function from the 
curvature perturbations. 


3.3 Curvature perturbations 

The spacetime primordial variations in the inhationary model have their origin in 
the quantum behavior of both the held and the spacetime itself. In section 3.1.2 
we studied the classical quantities (j) and in a background spacetime. Now we 
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want to study the quantum perturbations <50 and on the classical helds 

(j){t) (j){t) + 6(j){t, x), (3.25) 

x). (3.26) 

We only consider scalar perturbations since the tensor perturbation do not con¬ 
tribute to the bispectrum [41]. To begin with, we need to hnd the different scalar 
functions responsible for the general perturbations S(p and dg^^ty to the classical 
background. The hrst one is the scalar fluctuation to the scalar held S(p. To de¬ 
termine the other scalars from the metric perturbation we divide Sg^^iy into blocks 

Sgooit,x) 6goi{t,x) 

6gio{t,x) 6gii{t,x) 

We can now determine each of these scalar held. The easiest to determine is the 
one given by goQ which provides one scalar held. The term doi can be viewed as 
a spatial vector generated by taken the gradient of a scalar held B. The term 
5gij contains two scalar helds, one is its trace and the other one is given by ^ 

which contributes to 6gij by taking two spatial derivatives didj^. We thus have 

hve separate scalar held, and four of them correspond to the metric perturbations. 
However, in single held inhation only one of the hve scalar perturbations have a 
physical meaning [11]. 

In order to study the huctuations about the simple, spatially invariant inhationary 
background we follow the procedure described by Maldacena [17]. It is convenient 
to use the ADM formalism [17, 38] in which the metric (3.7) takes the form 

= —(A^^ — NiN^)dC -|- 2Nidtdx'^ + hijdx^dx^, 

where htj the spatial metric of the thre-dimensional hypersurfaces embedded in 
the full spacetime. N(t,x) is the lapse function and N’‘{t,x) the shift vector [26]. 
One advantage of this approach is that the helds N and W are both Lagrange 
multipliers with no underlying dynamics and whose equations of motion are easy 
to solve. Also these equations produce two constraints that reduce the number 
of independent scalar degrees of freedom by two [11, 13]. Notice that the spa¬ 
tial indices are contracted using the metric hij, for instance N'' = If^Nj. The 
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components of the inverse metric are 

= = = -N^n^). (3.29) 

Inserting the metric g^^i, into the full action of the theory Eq. (3.7) we have [13] 

S = ^J dtdx^VhN [i?® + N-^{EijE^^ - E^) (3.30) 

+N-^{^ - N^di(j)f - - 2I/(0)j , 

where is the 3-D Ricci scalar curvature constructed form hij and we have set 
Mp = 1 to simplify the calculations. We have introduced the spatial tensor Eij 
which is closely related to the extrinsic curvature Kij [26]. The latter determines 
the rate of change of the spatial metric hij, as we move along the timelike geodesics 
orthogonal to the hypersurfaces embedded in the spacetime [26]. The former is 
dehned by 

= I {hij - , (3.31) 

= ^ {hij - hikdjN^ - hjkdiN^ - N’^dkhij'^ , 

where is calculated with the curvature 3-D metric. The dehnition of E is 

E = RjE^h (3.32) 

To proceed with the ADM formalism we now vary the action with respect to the 
helds N and A^h Varying with respect to N we have 

fl{3) _ N-\EijE^^ - E^) - V-2(0 - N^di(j)f - EWi(j)dj(j) - 2V = 0. (3.33) 

From Eq. (3.10) we can replace the pontential for which we have not made any 
particular consideration, yielding 

^(3) ^ N-2^EijE^^ - E^) - - N^di(j)f - E^di^dj^ - 6H^ + (^^ = 0. (3.34) 

And with respect to W 


Vf [N-\Eh - SdE)] - N-\^ - NWj(j))di(j) = 0. 


(3.35) 
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A general parameterization of the held’s fluctuation is given by 

(j) = (l){t) + 5(l){t, x). (3.36) 

And of the scalar perturbations to the metric can be written down as [7, 13] 

7V= l + <h(f,f), (3.37) 

= S^^djB{t,x), (3.38) 

hij = a? [(1 + 2C(t, x)) 5ij + didj^]. (3.39) 


We can remove two of the hve scalar perturbations by a time redehnition and a 
shifting of the spatial coordinates x* ^ x* + S^^djf{t,x). We use this freedom to 
choose the coordinates such that <50 = 0 and ^ = 0. Hence the parameterizations 
on the helds are 


= 1 + ^{t, x), 

(3.40) 

W = 6^^djB{t,x), 

(3.41) 

hij = 

(3.42) 

0 = Ht), 

(3.43) 


where now 0 has no spatial dependence and ( has been slightly redehned for useful 
purposes when deriving the second and third order action of the theory. In this 
gauge the constraint equations Eqs. (3.34) and (3.35) are written 

- E ‘^) - (1 - iV - 2)02 _ g ^2 ^ ('3 44 ^ 

Vf [N-\E,^ -5,^E)] =0. (3.45) 


We now solve the constraint equations in terms of the scalar held huctuations. In 
these coordinates the scalar curvature is written 


i?(3) = -4a-^e-^< 


d,dX - -dXdX 


= -4a-^d,dX + 0((X, 


(3.46) 


where we have expanded ( up to hrst order. Similarly 


Eij = 


H + C- dXd^B 


didjB 


= a 


77(1 + 2C)5*, + C5,j - d,d,B 


(3.47) 

(3.48) 
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then to hrst order in ( 

EijE^^ -E^ = -QH^ - 12HC + AHdidjB + ... (3.49) 


In this way the constraint equation Eq.(3.44) can be written to hrst order as 


- 3H 


2H^-C 


- did^ [HBa-^C] + 0 ^ 4 ) = 0 . 


(3.50) 


And Eq. (3.45) takes the form 


29, 


2if4>-C 


0 . 


(3.51) 


Thus we have that 

^ = ~ + constant = (3-52) 

for which the constant is taken as zero in order to recover the original background 
metric if the perturbations are removed [13]. Inserting 4> into Eq. (3.50) we have 


9,9* [HB + a-^C] + = 0. (3.53) 

And solving for B 

with 

X = (3.55) 

where 9“^ is the inverse Laplacian and we have used the expression for e, Eq. 
(3.15). The dehnition of y given here differs from that of [13], but agrees with 
the one in Refs. [11, 41]. We have now solved the constraints N and iV* and have 
obtained expressions for them in terms of only one scalar held perturbation ( 


N 


TV* 


hij 


1 + 
9 


U 

2H' 

1 , 1 ■ 
—iijC + —X 




(3.56) 

(3.57) 

(3.58) 
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To obtain the two and three-point correlation functions we write the full action in 
the chosen gauge where 6(f) = 0 


S = 


1 

2 


dtdx^VhN 


rC) + 


E^) - 6H^ + {N + 


(3.59) 


Then we express all terms h, Etj ,... in terms of the physical scalar held (f 
and expand in terms of ( up to the desired order. 


Expanding to second order in the small huctuation ( we have 


^2 = Ml 


dtd^x 


- ae{dCf 


(3.60) 


where we have restore the Planck mass Mp. Expanding to third order 


S, = Ml 


dtd^x 


a?e 


aR CC + ae ((^C) - 2 aeC( 5 C)(ax) + —r]C ( (3.61) 
+ ^{dC){dx)d\ + + /(O^ 


where 


SL 

/(C) 


2a + Hd^x - , 

Tj 

—C + terms with derivatives on C 

4 


(3.62) 

(3.63) 


Here 6L/6(\i is the variation of the quadratic action with respect to (. In the next 
sections we will use the second and third order action expressions to calculate the 
power spectrum and bispectrum of perturbations, respectively. 


3.4 Power spectrum and Bispectrum 

Power spectrum For a long time, observations of huctuations in the CMB 
and LSS have been focused mainly on the Gaussian contribution as measured by 
the two-point correlations of the density huctuations or by its Fourier transform 
[ 10 , 11 ] 

= {27if‘^P^{k)6^^\ki + k2), 


(3.64) 
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where P({k) is the power spectrum of curvature perturbations dehned by 

Pdk) ^ (3.65) 

where Tg is the time at the end of inflation [40, 42], It has been precisely the 
study of the power spectrum which has settled the numerical value of many of the 
cosmological parameters in the ACDM model [1-3]. 


Bispectrum The three-point correlation function correlates density or tem¬ 
perature fluctuations at three points in space. Its Fourier transform correlates 
fluctuations with three wave vectors ki {i = 1,2,3) and is called the bispectrum 
given by the relation [10, 43] 

(^C{ki,t)C{k2,t)C{k3,t)'j = {27rfB(;{ki,k2,k3)d^^\ki + k2 + k3), (3.66) 

and it vanishes if the curvature perturbation is Gaussian [11, 40]. Therefore, 
deviations from Gaussianity can be determined by measuring A convenient 
way to do this is to parametrize the level of non-Gaussianity by using a Non- 
Linearity parameter /jvl [11, 43]. For instance the local model [44, 45] is dehned 
in terms of a local quadratic correction to the curvature perturbations given by 
[46] 

$(T) = ^{x)l + fNL [Hx)l - (<l>(f)L)'] , (3.67) 

where $(x)i is the linear gaussian part of the perturbation and characterizes 
the amount of non-Gaussianities. The non-linear coupling constant Jnl is also 
known as the bispectrum amplitude [43]. 

Note that the bispectrum is a function of the scales ki,k 2 , and k^, but since one 
is hxed by the function it only depends on two of them. These vectors form a 
triangle in Fourier space [1, 11, 13]. The dependence of on the scales is divided 
in two parts. The hrst one is called the shape of the bispectrum and corresponds 
to the dependence of 5^ on the momenta ratio k 2 /ki and k^/ki having the total 
momentum ktot = ki + k 2 + k^ hxed. The second one is called the running of 
the bispectrum corresponding to the dependence of Bc^ on the overall momentum 
ktot with hxed momenta ratio k 2 /ki and /cs/Zci. The most studied shapes are 
the local form or squeezed form [11, 41, 43], corresponding to triangles such that 
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fci -C ^2 = ^3 and the equilateral form corresponding to ki = k 2 = k^- In this 
thesis we will study both shapes thoroughly. 

The Planck constraints on primordial non-Gaussianity on the bispectrum ampli¬ 
tude for the local shape is [43] 


4Y' = 2.7±5.8, (3.68) 

and on the equilateral shape 

f^fateral ^ _^2 ± 75 . ( 3 . 69 ) 


Following the procedure in Ref. [41], where there is a held redehnition and some 
boundary terms are kept, the third order action Eq. (3.61) can be simplihed as 


S's D Mp / dt(fx 


-a^er/CC^ - -aer/C^^C 


(3.70) 


In the in-in formalism the above action can be interpreted in the interaction pic¬ 
ture. The Hamiltonian after changing to conformal time gives 


H,nt{r) = Ml 


d^x erja 




(3.71) 


The 3-point correlation function at a time Te after horizon crossing is then given 
by [11, 17] 


CL 


C(re, /Cl)C(re,/C2)C(R,^3) 


Cl) = —i 


C(r, fci)C(R, /s2)C(r, h). Hint 


(3.72) 

where according to the standard procedure Hint evolves the free theory vacuum 
state |0) to the interaction vacuum |r2) up to the time where the correlation func¬ 
tion is to be evaluated. After substitution of Eq. 3.71 the expression for the 
bispectrum [11, 47] is 


B^{kuk2,k3) = 2{2 tiYM\ 




C(a, fci)C(A, fc2)C(A,fc3) / dT7]ea\*{T,h) 


(3.73) 


X (2C'*(a k2)C{r, fcs) - klC{r, k2)C{r, k^)) 


+ the other two permutations of ki,k 2 , and fca. 
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where O' is the imaginary part and Te is some time sufficiently after Hubble cross¬ 
ing horizon [40, 42, 48] when the modes are frozen. To study the squeezed and 
equilateral limits of the bispectrum it is useful to use the dehnition function 
[49] given by 


FNiih, /C2, fcs; k*) 


10 {kik2k^)'^ 1 

3(27r)^ kl + kl + k^ P^{k^) 


B(^{ki, /C2, /cs), 


(3.74) 


where /c* denotes a pivot scale for the power spectrum [10, 41]. 




Chapter 4 


Effect of features in the potential 
of the inflaton 


In this chapter we study the effects on the power spectrum and bispectrum of 
primordial curvature perturbations produced by a discontinuity in the n — th order 
derivative in the inflaton potential and the suppression of the value of the inflaton 
held during inhation. 


4.1 Scale Dependent Features 

In the last decades, the information from the observations of huctuations in the 
CMB and the Large-Scale Structures have been extracted from the Gaussian con¬ 
tribution [1, 12, 17]. But any information contained in the departure from a perfect 
Gaussian held is not encoded in the power spectrum, thus it has to be extracted 
from measurements of higher-order correlation functions. From this measurements 
we can diherentiate inhationary models that can lead to very similar predictions 
for the power spectrum of primordial perturbations. Thus non-Gaussianity is a 
sensitive probe of the interaction of the helds driving inhation and therefore con¬ 
tains important information about the fundamental physics during inhation. 

At the end of section 3.2 we enumerate various mechanism to produce a large 
amount of non-Gaussianities. In this thesis we study a model where the slow- 
roll conditions are violated by the introduction of a “feature” in the potential or 
more generally in the Lagrangian of the inhaton held. In this way we may study 
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how the non-Gaussianities can discriminate between different inflationary models 
which can predict the same power spectrum. 

There are several other reasons that motivate the introduction of features, which 
have a long history [50] in the study of primordial fluctuations. For instance, it 
has been shown that these models can provide better hts to the apparent low 
multipole glitch at / ~ 20 — 40 in the study of the angular spectrum of the CMB 
radiation [51, 52]. Moreover, in the Planck 2013 results a new feature was found 
in the CMB spectrum at / ~ 1800 [10]. According to the Planck Collaboration, 
the features in the temperature power spectrum, particularly the broad dip at 
I ~ 1800, cannot be explained by the standard ACDM model. These departures 
from the best £t ACDM model spectrum could also be due to unknown system¬ 
atic effects into the hnal power spectra [10]. Only when the whole data from the 
mission can be analyzed this possibility will be investigated. To study the devi¬ 
ations of the primordial power spectrum from a smooth function various models 
have been proposed (See Ref. [53] for a thorough list and explanation of many 
different inflationary models). Some models can add a global oscillation, a local¬ 
ized oscillation, or a cutoff to the large scale spectrum, namely the Wiggles, the 
step-inflation, and the Cutoff models, respectively [10]. For instance, in the case of 
step-inflation it was shown in Ref. [54] that a large number of “sudden downward” 
steps in the inflaton potential can be obtained from a class of models derived from 
supergravity theories. In these models the symmetry breaking phase transition of 
a held coupled to the inhaton “naturally” gives rise to the steps, since the mass 
changes suddenly when each transition occurs. 

In this thesis we are interested in the effects that an n-th order discontinuity can 
produce (on the power spectrum and bispectrum) [10, 40, 42, 50, 55] and on the 
consequences of having a term which suppresses the value of the scalar held during 
inhation. This corresponds to a non-Iocalized feature in which the ehect of the 
suppression is higher at the end of inhation. 
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4.2 The Model 


We consider a single scalar field 0 coupled to a potential given by 


v{4>) 


Vob + (j) > 00 , 

Voa + \nF(tF + AA0 0 < 00, 




(4.1) 


A0 = 0”e \ ^0 ' , 


where m is the mass infiaton. Vofe and Voq are the vacuum energies before and 
after the transition, respectively. The transition time is Tq and we define 0(ro) as 
00 . In Eq. (4.1) n is a positive real number and A is a model parameter whose 


units in terms of Planck mass are We will see below that this feature in 


the potential induces an oscillatory ringing in the power spectrum of curvature 
fluctuations [55]. The case for n = 2 corresponds to a sharp feature in the infiaton 
mass [41] as we will see below. 

The equations of motion (3.9) and (3.10) in terms of conformal time r (2.18) take 
the form 



(4.2) 


4>" + 2-4>' + a'‘dsV = 0, 


(4.3) 


a 


where primes indicate derivatives with respect to conformal time. Imposing con¬ 
tinuity conditions on the potential at 0o implies that Voa = Vob — A0o. We assume 
that this change in the potential energy is small [40, 55]. In the next section we 
will find an analytic solution for a and 0 in terms of r. 

4.2.1 Analytic solutions to the equations of motion 

In order to obtain analytic solutions for the equations of motion we assume that 
V (0) is dominated by the vacuum energy Vq. In this case H is constant and the 
scale factor may be approximated by that of a pure de Sitter space 



(4.4) 
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We now proceed to obtain the analytic solution of the held. Before the transition 
Eq. (4.3) is written as 

(4.6) 


+ 2—(f) + a m (j) = 0, 
a 


which has the solution 


Mn = 4>t<^oF + <t>t «(!■)*, 


(4.6) 


where the slow-roll regime corresponds to 0^ =0 and is dehned as 




2 m 


(4.7) 


After the transition Eq. (4.3) is written as 

a' 

(j)" + 2—0' + a? (rn^cj) + = 0, 

Qj 

where 

d4,t\(p = + 24," 


(4.8) 


(4.9) 


Eq. (4.8) does not have an exact analytic solution. To obtain an analytic approx¬ 
imation we expand 0 to second order around the transition 


Hr) = 00 + 0'(ro)(r - Tq) + ^0"(ro)(r - Tq)^ + C>(3), 
and then insert this expansion into 0</,A0 to obtain 


(4.10) 


+ 2—4' + 

a 


m (j) A0o" ra0o 


(4.11) 


(n+ l)(n - 2)4'(to)(t - To) + -(n + l)(n - 2)4"(to)(t - To)^ | = 0, 


where we have neglected the term ?7,(?7, — 4)0'(ro)^ with respect to (n — 2)0(ro)0"(ro) 
since 


n(n — 4)0'(ro)^ n(n — 4) m^ 


in 


< 1 , 


n ^ 2 


(4.12) 


2)0(ro)0"(ro) 3{n - 2) 

in the last step we use the fact that m H otherwise the scalar held oscillates 
with a frequency proportional to its mass^ then it would not correspond to a 


^This can be seen in Eqs. (4.6) and (4.7). See also Ref. [7] exercise 6.7. 
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slow-roll regime. The solution is thus given by 


0a('r) = + 0a’ (r - To) + 0a’ (t - Tof + 0+a(r)^'^ + 0„ a(r)^ , (4.13) 


where 




-A 0 O 


n—2 


[n{m^ — 2iJ^)0o 


+2(n + l)(u - 2)hfVoro - (u + l)(u - 2)ifVoVo'] , 

— Xrh 

+!)(«-2)4. 


(4.14) 


(4.15) 

(4.16) 


The quantities 0a’ {i = 1, 2, 3.) are in terms of 0o, 0o = 06(^0), and 0o = 0f(('ro) for 
which we already have a solution. The constants of integration 0^ are determined 
by imposing the continuity conditions on 0 and 0' across Tq 


^>a 


±a(ro) 


-A± 


(A- - A+) 


A^0O + 00 T) + 

{n + l)(n — 2) 
(m 2 - 2H‘^) 


A 0 O 


n—2 




nA^ 0 o 


(4.17) 


((m= + 2H^-Xf )4>'aT„ - 


Thus our complete analytic approximate solution is 

0(r)analytic = (j)b{T) + 0a('r)6'(0O “ 0), (4-18) 

where 0b and 0a are dehned in Eqs. (4.6) and (4.13), and 6*(0) denotes the Heav¬ 
iside step function. 


We choose the model parameters as [41] 

m = 6 X IQ-^Mpi, H = 2x lO’^Mp^, 0+ = lOMpp (4.19) 

From now on we set Mpi = 1, thus all quantities will be given in Planck units; for 
instance the transition time is tq = —Mpl = —1. 
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{n=3, A=2.4xl0-‘*) 

0 



Figure 4.1: Evolution of (p for n = 3 and A = 2.4 x 10“^® in terms of confor¬ 
mal time. The blue and dashed-red lines represent the numerical and analytic 
results, respectively. The agreement between these two solutions is observed. 
All quantities are in Planck units. 

4.2.2 Background solution and slow-roll parameters 


In this section we show that our analytic solution is a good approximation for 
the held inhaton under the de Sitter approximation. To do this we compare the 
numerical solution for a and 0 [Eqs. (4.2) and (4.3)] with the analytic solution 
Eq. (4.18) in the de Sitter approximation. In the following plots all quantities are 
in Planck units. It is important to remember that the slow-roll parameters are 
dimensionless. 


In Fig. 4.1 we plot the numerical and analytical evolution of 0 as a function 
of conformal time. The evolution corresponds to a slow-roll regime. There is 
an agreement between the numerical and analytic solutions for the scalar held. 
Fig. 4.2 shows that the percent error between these two solutions is very small: 
after 60e-folds this error is only of 0.4%. We also found agreement between the 
numerical and analytic solution for other values of n and A. We can conclude then 
that Eq. (4.18) is a good approximation to the background solution, thus in the 
following sections we use the analytic solution Eq. (4.18) to obtain the analytic 
results for the plots. 


The slow-roll parameters Eqs. 
the form 

e 


(3.13) and (3.14) in terms of conformal time take 


H' _ e' 


(4.20) 
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{n=3, fc2.4xl0-‘*) 

ioo%|— - - 



Figure 4.2: Plots of percent error between the numerical and analytic solutions 
for the field inflaton (j) for n = 3 and A = 2.4 x 10“^®. Although this percent 
error is increasing, after 60e-folds the error is only of 0.4%. 


We approximate the analytic solutions for the slow-roll parameters after the tran¬ 
sition as 


(^) 



(A+0Xr)^’" + A 0„a(r)^ ) 


2 


(4.21) 


^a(r) 


2A- 


+0-a(r)^' 


(4.22) 


For comparison, before the transition p = 2A’''. Due to the step in the potential 
the inflaton suffers a strong instant acceleration [55]. But still for realistic models, 
it is required that the step be approximately less than 1% of the overall height of 
the potential [42, 55]. Thus, although e grows after the transition, we still have 
that at all times cfF D(0) and e <C 1. Furthermore, the step causes a large 
change in V” and from Eq. (3.20) we have that rj oc V"{(j)) then t] will suffer a 
dramatic change [12, 55]. This is the primary source of large non-Gaussianities 
in the discontinuity of the 2-nd order derivative of the potential as can be seen 
in Eq. (3.73). Figs. 4.3 and 4.4 shows the evolution of e and rj for n = 3 and 
A = 2.4 X 10“^® over the transition. The analytic Eqs. (4.21) and (4.22) are good 
approximations for the slow-roll parameters. Similar results are obtained for other 
values of n and A in which there is also an agreement between the numerical and 
analytic solutions. 
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{n=3, A=2.4xl0“‘*) 



Figure 4.3: Plots of e in terms of conformal time showing the evolution over 
the transition for n = 3 and A = 2.4 x 10“^®. The blue and dashed-red lines 
represent the numerical and analytic results, respectively. 


{n=3, A=2.4xl0-‘*) 

f] 



Figure 4.4: Plots of the slow-roll parameter in terms of conformal time for 
n = 3 and A = 2.4 x 10“^®. This violation of 1 is the primary source of 
the large non-Gaussianities in the discontinuity of the 2-nd order derivative of 
the potential, as the leading term in the Fjvl function Eq. (3.74) is of order 
rj. The blue and dashed-red lines represent the numerical and analytic results, 
respectively. 

4.3 Perturbations 


From the action Eq. (3.60) we obtain the linear equation of motion for the curva¬ 
ture perturbation ( as 


dt 





d'X 

dx'-dx^ 


0 . 


(4.23) 


In Fourier space and in terms of conformal time the previous equation is written 


in the form 



0 , 


(4.24) 
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where z = a\ple and k is the comoving wavennmber. Towards a better under¬ 
standing on the evolution of the curvature perturbation it is convenient to dehne 
a gauge invariant potential [41, 55] 


Mfc(r) = 2;(r)C(r, k) 


(4.25) 


Then Eq. (4.24) takes the form^ 


tlh 




Mfc 0, 


(4.26) 


for which the dependence on time has not been written explicitly. Let us now 
consider the limits ^ z”/z and /c^ <C z" jz. In the hrst limit 


Uk oc 


(4.27) 


which is the free held solution. In the second limit, we have 


Uk oc z. 


(4.28) 


which means that the curvature perturbation ( is constant 

Uk 

( = — oc constant. (4.29) 

In the slow-roll regime the term z” j z can be written approximately as z”j z = 
2a^H^ . The physical wavelength of the perturbation is proportional to a/k and 
the physical horizon during inhation is given by H~^. We thus have that in the 
hrst limit the modes are well within the horizon 

k 3> aH modes inside the horizon. (4.30) 

While in the second limit 

k aH modes outside the horizon. (4-31) 

The transition time tq sets the scale of transition at ko = —1/to = 1. At the time 
of transition there is a discontinuity in 0" which implies that z” contains a Dirac 
delta function [56]. To evaluate this discontinuity we calculate the contribution of 

^The definition of Uk in Ref. [55] differs from our definition since in our case 2 is a positive 
function. 
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the delta function in z" as 


Do = lim / —dr = 


5^0 


TQ—S 


Z 


^ bao - = -nXa{To) 

ro 


in—l 


2^0 


(4.32) 


Our choice of vacuum implies that the initial conditions for the mode function 
and its derivative are given by the Bunch-Davies vacuum at early times r —?■ cxd 
[17, 41] 

n(r,/c) 1 \ 


ri 7^ ^ (v{T,k)\ 

C('r, k) = and C (b k) = dr [ , . 

z{t) \ z{t) J 


where 


n(r, k) = 


~,—ikT 




kr 


(4.33) 

(4.34) 


An analytic solution to the mode function is given by [41] 


C(r, k) = 


v{T,k) 

y/M-r) 

v{T,k) 

\/2eCfc) 

v{T,k) 


Mpair) 


k < ko and r < Tk 
k < ko and t > 
k > ko and r < tq 


\/2e(T) 

a(k)v{T,k)+f3{k)v* (T,k) ; ^ ; j ^ ^ 

k > ko and tq < r < Tk 
c.{k)v(r,k)+mv*{Tk) > fco and r > Tfc, 


where we have restored the Planck mass and 


(4.35) 


a{k) = 1 + iDo\v{To,k)y and /3(/c) = —iDov{To,kY, 


(4.36) 


and Tfc = —l//c is the horizon crossing time for that mode. We will see in Sec. 4.4 
that this solution is a good approximation for the power spectrum. 


4.4 Results and Discussions 

4.4.1 Background solution and slow-roll parameters 

The numerical results are presented in Figs. 4.5-4.7. The time evolution of the 
inflaton is plotted in Fig.4.5 for different types of features of the potential. These 
effects on the slow-roll parameters can also be seen in the Fig. (4. 6-4. 7). When A 
is kept constant, larger values of n tend to produce larger increments of both e 
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Figure 4.5: On the left the evolution of (f> is plotted for A = 3.8 x 10“^® 
and n = 2/3(blue lines), n = 3 (red lines), and n = 4 (green lines). On the 
right the evolution of (p is plotted for n = 3 and A = 6.0 x 10“^® (blue lines), 
A = 1.2 X 10“^® (red lines), and A = 2.4 x 10“^® (green lines). 


E E 




Figure 4.6: On the left the slow-roll parameters e is plotted for A = 3.8 x 10“^® 
and n = 2/3 (blue lines), n = 3 (red lines), and n = 4 (green lines). On the 
right the slow-roll parameter e is plotted for n = 3 and A = 6.0 x 10“^® (blue 
lines), A = 1.2 X 10“^® (red lines), and A = 2.4 x 10“^® (green lines). 


and rj, while when n is kept constant, larger values of A tend to produce larger 
increments of both e and rj. 


4.4.2 Power spectrum 


In section 3.4 we saw that the power spectrum of curvature perturbations was 
dehned as Eq. (3.65) 

9P 

(4.37) 


0^3 


(2ir) 

For small scales, k ':$> ko, there is a simple expression for the power spectrum [56] 

1 1/2 


p - G/2 ^ a{Tk)H{Tkf 
27r|0'(rfe)| 


1 + ^ sin(2fcro) + ^ [^ + cos(2/cro)] 


(4.38) 
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Figure 4.7: On the left the slow-roll parameter rj is plotted for A = 3.8 x 10“^® 
and n = 2j3 (blue lines), n = 3 (red lines), and n = 4 (green lines). On the right 
the slow-roll parameter rj is plotted for n = 3 and A = 6.0 x 10“^® (blue lines), 
A = 1.2 X 10“^® (red lines), and A = 2.4 x 10“^® (green lines). This violation of 
the slow-roll condition is the primary source of the large non-Gaussianities in 
our model since the leading term in the Fjvl function is of order r]. 





Figure 4.8: Plots considering two different analytic expressions for the power 
spectrum: the first one (dashed-red) is given by inserting the analytic approxi¬ 
mation (4.35) into the definition of Eq. (4.37) and the second one (dotted- 
green) is given by Eq. (4.38). These analytic expressions are compared with the 
numerical result (blue lines). The parameters used in the three plots are n = 3 
and A = 2.4 X IQ-^®. 


In this section we use this expression along with Eq. (4.37) to compare the power 
spectrum with the numerical result. In Fig. 4.8 we compare two different ana¬ 
lytic expressions for the power spectrum. In the hrst one we insert the analytic 
approximation (4.35) into the dehnition of the power spectrum Eq. (3.65). For 
the second one we use the analytic expression (4.38). To plot these functions we 
use the analytic solution of the scalar held Eq. (4.18). The analytic expressions 
are compared with the numerical result that we obtained from our code. The 
parameters used in the plots are n = 3 and A = 2.4 x 10“^®. As it can be seen 
in Fig. (4.9) small scales modes which are sub-horizon at the time of when the 
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Figure 4.9: The numerically computed |Cfc| is plotted as a function of the 
number of e-folds N after the time of the feature . The left plots are for A = 
3.8 X 10“^® and n = 2/3 (blue lines), n = 3 (red lines), and n = 4 (green lines). 
The right plots are for n = 3 and A = 6.0 x 10“^® (blue lines), A = 1.2 x 10“^® 
(red lines), and A = 2.4 x 10“^® (green lines). All the plots are for the short 
scale mode with k = IOO/cq which is sub-horizon when the feature occurs. 





Figure 4.10: The power spectrum of primordial curvature perturbations 
plotted for different types of features. For the left plots A is constant, A = 

3.8 X 10“^®, and n = 2/3 (blue lines), n = 3 (red lines), and n = 4 (green 
lines). For the right plots n is constant, n = 3, and A = 6.0 x 10“^® (blue lines), 

A = 1.2 X 10“^® (red lines), and A = 2.4 x 10“^® (green lines). 

feature occurs are affected by the feature, while modes that had already left the 
horizon are unaffected. In £g.(4.10) the power spectrum of primordial curvature 
perturbations is plotted for different types of features. 


4.4.2.1 Models with the same spectrum 

We have found an interesting case in which the the power spectrum for different 
values of n and A can be the same. For instance for (ni, Ai) and ( 77 , 2 , A 2 ) we can 
have 


Pc^ik, rii, Ai) = P^(/c, 712, A 2 ), 


(4.39) 
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Figure 4 . 11 : Numerical plots of for three different values of n and A for 
which the power spectrum is exactly the same. The parameters used are (blue) 
re = 2/3 and A = 2.3 x 10 “^^, (dashed-red) re = 3 and A = 2.4 x 10 “^®, and 
(green) re = 4 and A = 1.8 x 10 “^®. 



for all values of k. Those values correspond to setting 

nA(/>Q~^ = constant. (4.40) 

In this case Dq is the same for all models with the same set of parameters Eq. 
(4.19). The interesting fact about this is that the degeneracy can only be broken at 
the bispectrum level as we will see below. The results are shown in Fig. 4.11 where 
we use three different values for n and A and obtain the same power spectrum. 


4.4.3 Bispectrum 

The feature in the potential generates large gaussianities for the cases of large n 
or large A as is expected from Eq. (3.73) given by 


B^{ki,k2,k3) = 


3/1^2^ 


C,{re,ki)C,{Te,k2)C,{re,kY / drpeaXY'^^h 


'TO 


(4.41) 


X (2C"(t, i: 2 )C'*(T, is) - kjC{T,k2)C{r,h)) 

the other two permutations of ki,k 2 , and k^, 


since as we mentioned above as the inflaton crosses the feature r] increases dra¬ 
matically. We evaluate the integral from tq —)> Tg since for single held inhationary 
models ( —?> constant after Hubble crossing horizon as it freezes out, while for early 
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Figure 4.12: The squeezed limit of the numerically computed bispectrum 
FNL{ko/500, k, k) in plotted for a large scale ko/bOO. On the right we keep A 
constant, A = 3.8 x 10“^®, while n = 2j3 (blue lines), n = 3 (red lines), and 
n = 4 (green lines). On the left we keep n constant, n = 3, while A = 6.0 x 10“^® 
(blue lines), A = 1.2 x 10“^® (red lines), and A = 2.4 x 10“^® (green lines). 
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Figure 4.13: The squeezed limit of the numerically computed bispectrum 
FArL(A:, 1000/co, 1000/co) is plotted for a small scale IOOO/cq- On the left A is 
constant, A = 3.8 x 10“^®, while n = 2/3 (blue lines), n = 3 (red lines), and 
n = 4 (green lines). On the right n is constant, n = 3, while A = 6.0 x 10“^® 
(blue lines), A = 1.2 x 10“^® (red lines), and A = 2.4 x 10“^® (green lines). 


times C oscillates rapidly, so its contribution to the integral tends to 

cancel^ [40, 48]. In this way the integral is dominated by the perturbations leaving 
the horizon. The numerical results for different values of n and A are plotted in 
Figs. 4.12 - 4.15 for large and small scales in the squeezed and equilateral limits. 


4.4.3.1 Analytical approximation for the bispectrum 

In order to obtain an analytic approximation for the bispectrum we use the analytic 
expressions for the perturbation Eq. 4.35 and for the slow-roll parameters Eqs. 
4.21 and 4.22. We also use a series of approximations such as hxing e and rj constant 
after and before horizon crossing for those modes e(rfc) and ri{Tk), except in the 
time interval around when the feature occurs, when in fact their variation cannot 
be neglected. We also consider that A"*" <C A“ throughout the calculation. Another 


^See Refs. [40, 48] for a detailed treatment of the integral for times r < tq. 
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Figure 4.14: The equilateral limit of the numerically computed bispectrum 
FNL{k,k,k) in plotted for large scales. On the right we keep A constant, A = 
3.8 X 10“^®, while n = 2/3 (blue lines), n = 3 (red lines), and n = 4 (green 
lines). On the left we keep n constant, n = 3, while A = 6.0 x 10“^® (blue lines), 
A = 1.2 X 10“^® (red lines), and A = 2.4 x 10“^® (green lines). 
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Figure 4.15: The equilateral limit of the numerically computed bispectrum 
F/vl(/c, k, k) is plotted for small scales. On the left A is constant, A = 3.8x 10“^®, 
while n = 2/3 (blue lines), n = 3 (red lines), and n = 4 (green lines). On the 
right n is constant, n = 3, while A = 6.0 x 10“^® (blue lines), A = 1.2 x 10“^® 
(red lines), and A = 2.4 x 10“^® (green lines). 


approximation considers the fact that the modes freeze out after horizon crossing 
then after that the modes are constant and can be taken out of the integral. The 
analytical approximations for the bispectrum are shown in Figs. 4.27 - 4.23 and 
are in good agreement with the numerical results. 


Squeezed limit 

• In the squeezed limit and for small scales with ko ^ ki k 2 = = k we 

found an analytic expression for in terms of Dq given by 


F'; 




NL 


-4 


ko 


cos 


2k + ki 
kn 


2k + kiki f 2k + ki 

—I - 

ko ko 

Do 


ko 


■ cos 


ko 
2k — ki 
ko 


(4.42) 


where SL stands for squeezed limit and the sign > for small scales. 
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Figure 4.16: The squeezed limit of the numerically computed (blue lines) 
and analytic (red lines) bispectrum F;vl(A:o/500, k, k) is plotted for a large scale 
feo/500. Prom left to right we keep A constant, A = 3.8 x 10“^®, while n = 3, 
and n = 4, respectively. 


• In the squeezed limit and large scales with ki k 2 = ks = k ^ ko the 
analytic expression for all n is 

F-<W = -(n-2)|5(8+^)(A)\ (4.43) 

where the sign < means for large scales. In this case the sign of Fjvl is 
determined by the factor (n — 2). 


Equilateral limit 


• In case of the equilateral limit and small scales with ko ki = k 2 = ks = k 
the analytic expression is written as 


t7EL> 
^ NL 




-4 


h 


■ sm 


3fc\ 


3--^ 
3 ko 


““ (I) 


(4.44) 


where EL stands for equilateral limit. 


• In the equilateral limit and for large scales with ki = k 2 = ks = k ko the 
analytic expression is given by 


F^L^ik) 


-in-2) 


4/^0 

3ko 



(4.45) 


where again the sign of the Fnl function is determined by n. 
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Figure 4.17: The squeezed limit of the numerically computed (blue lines) 
and analytic (red lines) bispectrum F]\fL{ko/500, k, k) is plotted for a large scale 
fco/500. From left to right we keep n constant, n = 3, while A = 6.0 x 10“^®, 
A = 1.2 X 10“^®, and A = 2.4 x 10“^®, respectively. 
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Figure 4.18: The squeezed limit of the numerically computed (blue lines) and 
analytic (red lines) bispectrum FArL(A:, lOOOfeo, lOOOfco) is plotted for a small 
scale lOOOfeo. From left to right we keep A constant, A = 3.8 x 10“^®, while 
n = 3, and n = 4, respectively. 
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Figure 4.19: The squeezed limit of the numerically computed (blue lines) and 
analytic (red lines) bispectrum FAri(A:, IOOO/cq, lOOOfco) is plotted for a small 
scale lOOO/co- From left to right we keep n constant, n = 3, while A = 6.0 x 10“^®, 
A = 1.2 X 10“^®, and A = 2.4 x 10“^®, respectively. 




Figure 4.20: The equilateral limit of the numerically computed (blue lines) 
and analytic (red lines) bispectrum F]\fL{k, k, k) is plotted for large scales. From 
left to right we keep A constant, A = 3.8 x 10“^®, while re = 3, and re = 4, 
respectively. 
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Figure 4.21: The equilateral limit of the numerically computed (blue lines) 
and analytic (red lines) bispectrum F]\fL{k, k, k) is plotted for large scales. From 
left to right we keep n constant, n = 3, while A = 6.0 x 10“^®, A = 1.2 x 10”^®, 
and A = 2.4 x 10“^®, respectively. 



Figure 4.22: The equilateral limit of the numerically computed (blue lines) 
and analytic (red lines) bispectrum F]\fL{k, k, k) is plotted for small scales. From 
left to right we keep A constant, A = 3.8 x 10“^®, while n = 3, and n = 4, 
respectively. 
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Figure 4.23: The equilateral limit of the numerically computed (blue lines) 
and analytic (red lines) bispectrum F]\fL{k, k, k) is plotted for small scales. From 
left to right we keep n constant, n = 3, while A = 6.0 x 10“^®, A = 1.2 x 10“^®, 
and A = 2.4 x respectively. 

4.4.3.2 Models with the same bispectrum 

We have also found some interesting effects in the F^l functions Eq. (3.74) when 
the appropriate choice of parameters give the same spectrum as shown in Fig. 
4.11. In this case in the squeezed and equilateral limits at large scales the Fnl 
functions are positive for n < 2 and negative for n > 2. While in both limits for 
small scales the Fnl functions could be the same when n and A fulhll Eq. 4.40. 
This corresponds to a break in the degeneracy since, as we saw above, the power 
spectrum was indistinguishable for different values of n and A. While in case of 
the squeezed and equilateral limits for small scales and for different values of n 
and A we found that the Fnl functions can be the same. 

In the squeezed limit and for small scales we found that Fnl is the same for 
different values of n and A as can be seen in Fig. 4.24. Fig. 4.25 shows the 
comparison between the numerical result and the analytic expression. For n > 2 
we have that Fnl < 0 as can be seen in Fig. 4.26 for which the numerical result 
is plotted for different values of n and A. The comparison between the numerical 
and analytic results in cases of n < 2 and n > 2 are shown in Figs. 4.27 and 4.28, 





















































































































































Chapter 4. Features in the potential 


54 


^NL 



Figure 4.24: Plots of Fjvl(A:, IOOO/cq, lOOOfco) in the squeezed limit of small 
scales and for three different values of n and A. The parameters used are (blue) 
n = 2l2> and A = 2.3 x 10“^^, (dashed-red) n = 3 and A = 2.4 x 10“^®, and 
(dashed-green) n = 4 and A = 1.8 x 10“^®. In this case the functions are 
the same for the appropriate values of n and A for which the bispectrum is the 
same. 
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Figure 4.25: Numerical (solid-line) and analytic (dashed-line) plots of 
FArL(fc, 1000/co, 1000/co) in the squeezed limit of small scales and for n = 
3 and A = 2.4 x 10“^®. This hgure shows the agreement between the numerical 
and analytic results. 


respectively. It can be seen that the analytic expression Eq. (4.43) determines the 
change of sign in Fnl correctly. 

In case of the equilateral limit and small scales we again found, Fig. 4.29, that 
the Ftvl functions are the same for different values of n and A. In the equilateral 
limit and for large scales the sign of the F/vl function is determined again by n. 
Figs. 4.30 and 4.31 show this change. The analytic expression is compared with 
the numerical results in Fig. 4.32. 
































Chapter 4. Features in the potential 


55 


^NL 



Figure 4.26: Numerical plots of FjYi(ko/500, k, k) in the squeezed limit of 
large scales and for three different values of n and A. The parameters used are 
(blue) n = IjZ and A = 2.3 x 10“^^, (red) n = 3 and A = 2.4 x 10“^®, and 
(black) re = 4 and A = 1.8 x 10“^®. In this case the functions are different 
for the particular values of re and A, although the power spectrum for this values 
is the same. 
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Figure 4.27: Numerical (solid-line) and analytic (dashed-line) plots of 
FNL{ko/5Qt)^ k, k) in the squeezed limit of large scales and for re = 2/3 and A = 
2.3 X 10“^^. This figure shows the agreement between the numerical and analytic 
results. 


The feature in the potential studied in this thesis can generate the same power 
spectrum for an appropriate choice of values on n and A. Although at the bis¬ 
pectrum level this degeneracy is removed. In this case we can see that at large 
scales and for values of re < 2 the Fj^l functions are positive in the squeezed and 
equilateral limits. While in the squeezed and equilateral limits for small scales and 
for different values of re and A the Fnl functions are the same. 

We summarized our results in Tables 5.1 and 5.2. For large scales we have k < ko 
while for small scales k > k^. The “Equal” indicates that the Fnl functions can 
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Figure 4.28: Numerical (solid-line) and analytic (dashed-line) plots of 
FNL{ko/500, k, k) in the squeezed limit of large scales and for n = 3 and A = 
2.4X10“^®. This hgure shows the agreement between the numerical and analytic 
results. 
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Figure 4.29: Plots of Fj\fL{k,k,k) in the equilateral limit of small scales and 
for three different values of n and A for which the bispectrum is exactly the 
same. The parameters used are (blue) n = 2/3 and A = 2.3 x 10“^^, (dashed- 
red) n = 3 and A = 2.4 x 10“^®, and (dashed-green) n = 4 and A = 1.8 x 10“^®. 
The errors marks in the plots correspond to a small path in the interpolation 
function to be plotted. In this case the F/vl functions are the same for the 
particular values of n and A. 


be the same for an appropriate choice of the values of the parameters for which 
the power spectrum is also the same. 
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Figure 4.30: Numerical (solid-line) and analytic (dashed-line) plots of 
Pwrik, k, k) in the equilateral limit of large scales and for n = 2/3 and A = 
2.3 X 10“^®. This hgure shows the agreement between the numerical and ana¬ 
lytic results. 



Figure 4.31: Numerical (solid-line) and analytic (dashed-line) plots of 
FNL{k,k,k) in the equilateral limit of large scales and for n = 3 and A = 
2.4 X10“^®. This hgure shows the agreement between the numerical and analytic 
results. 



k < ko 

k > ko 

Pc 

Equal 

Equal 


Table 4.1: Equality/Inequality of the power spectrum for an appropriate 
choice of n and A for which the spectrum can be the same. 



SL fc < fco 

SL k > ko 

EL k < ko 

EL k > ko 

Fnl 

Not equal 

Equal 

Not equal 

Equal 


Table 4.2: Equality/Inequality of the F/vl functions for an appropriate choice 
of n and A for which the spectrum can be the same. SL stands for squeezed 
limit and EQ for equilateral limit. For large scales we have k < ko while for 
small scales k > kQ. 





















Chapter 4. Features in the potential 


^NL 



Figure 4.32: Numerical plots of Fj^L{k,k,k) in the equilateral limit of large 
scales and for three different values of n and A. The parameters used are (blue) 
n = 2/3 and A = 2.3 x 10“^®, (red) n = 3 and A = 2.4 x 10“^®, and (black) 
n = 4 and A = 1.8 x 10“^®. In this case the functions are different an 
appropriate choice of n and A, although the power spectrum for this values is 
the same 









Chapter 5 


Conclusions 


In this thesis we studied the generation of primordial perturbations in the con¬ 
text of inflationary models with features in the inflaton potential. In Chapter 1 
we introduced a justihcation for the study of non-Gaussianities of cosmological 
perturbations. In Chapter 2, we reviewed the theoretical framework, assumptions 
and observational evidence which contributed to the establishment of the Big Bang 
model and highlighted its successes and problems. We then moved on to an ex¬ 
tension of the Big Bang model in order to account for solutions to the flatness and 
horizon problems as well as the formation of primordial perturbations. The result¬ 
ing model is presented in Chapter 3 as the Inflationary model, in which a scalar 
field minimally coupled to gravity drives an accelerated expansion of the universe 
solving the Big Bang problems and giving rise to the seeds of primordial perturba¬ 
tions. In this chapter we also introduced the calculation of the n-point correlation 
functions which relate the fluctuations of the spacetime metric and the inflaton at 
different places of the early universe. As explained in Section 3.2 many single-field 
slow-roll models of inflation predict small amounts of non-Gaussianity (correla¬ 
tion functions of order larger than 2), even though there is still a mechanism for 
producing a large and detectable amount of non-Gaussianity, which corresponds 
among others, to the violation the slow-roll conditions [1]. In Section 3.3 we re¬ 
viewed the Maldacena’s calculation [17] of the second and third order terms in 
the inflationary action. The power spectrum and bispectrum of the primordial 
curvature perturbations are defined in Section 3.4. 

In Chapter 4 we studied the effects on the power spectrum and bispectrum of 
primordial curvature perturbations produced by a discontinuity in the n-th order 
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k < ko 

k > ko 


Equal 

Equal 


Table 5.1: Equality of the power spectrum for an appropriate choice of n and 
A for which the spectrum can be the same. 


derivative in the potential of the inflaton. In Section 4.1 we reviewed the im¬ 
portance and justihcation of these features on the potential with respect to the 
observational data. As we mentioned there these models can provide better hts to 
the apparent low multipole glitch at / ~ 20 — 40 in the study of the angular spec¬ 
trum of the CMB radiation [51, 52]. One important characteristic of our model is 
that it can generate large non-Gaussianities through the violation of the slow-roll 
conditions. We also found that each different type of feature has distinctive effects 
on the spectrum and bispectrum of curvature perturbations which depend both 
on the order n of the discontinuous derivative and on the amplitude A of disconti¬ 
nuity. We have found that the spectra of primordial curvature perturbations have 
an oscillatory behavior for small scales. Also in the squeezed and equilateral limits 
for small scales the bispectrum has an oscillatory behavior whose phase depends 
on the parameters determining the discontinuity, and whose amplitude is inversely 
proportional to the scale. 

One of the most important results of our work is that the features in the potential 
can also generate the same power spectrum for an appropriate choice of values on 
n and A Eq. 4.40, although at the bispectrum level this degeneracy is removed. In 
this case in the squeezed and equilateral limits at large scales the Fnl functions 
are positive for n < 2 and negative for n > 2. While in both limits for small 
scales the Fnl functions could be the same when n and A fulhll Eq. 4.40. The 
importance of this result is that if we would only considered the power spectrum 
of curvature perturbations different models of inflation could be indistinguishable. 
It is then important to study the bispectrum of primordial perturbation in order 
to distinguish between several models of inflation [10]. We show again our results 
in Tables 5.1 and 5.2. For large scales we have k < ko while for small scales k > k^. 
The “Equal” indicates that the Fnl functions can be the same. 

Our numerical results of the power spectra are in agreement with the observational 
data from WMAP and Planck [10, 57]. The next step in the study of our model is 
to contrast the numerical results of the bispectra with the observational data. The 
interpretation of the CMB polarization from Planck data [10] or other experiments 
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SL fc < fco 

SL k > ko 

EL k < ko 

EL k > ko 

Fnl 

Not equal 

Equal 

Not equal 

Equal 


Table 5.2: Equality/Inequality of the T/vl functions for an appropriate choice 
of n and A for which the spectrum can be the same. SL stands for squeezed 
limit and EQ for equilateral limit. Eor large scales we have k < ko while for 
small scales k > kQ. 


such as Bicep2 have not been released yet. Even thongh onr results of the bispectra 
agreed with the previous results of WMAP [57] and the preliminary results from 
Planck [10, 43]. After comparing the data we could discriminate between different 
inflationary cosmological models in order to yield models that better describe the 
evolution of our universe. 

Also up to this point our analysis has been focused on the study of the two- 
and three-point correlation functions of primordial perturbations. We have been 
able to characterize different inflationary model according to their predictions 
and comparison with the observational data present up-to-date. But there is 
still a long way to go. Since different models of inflation can produce the same 
spectra and bispectra (Figs. 4.24 and 4.29) we should move on to the study 
of higher order correlation functions such as the four-point correlation function 
known as the trispectrum. Then we could distinguish between different models 
of inflation. To achieve this task properly our numerical code would have to be 
improved. Moreover the resolution of experiments which detect tiny correlations 
of temperature at different points in the sky would also have to be improved if we 
really want to compare the observational data with the predictions of inflationary 
models which produce non-Gaussianities. 


It is an exciting time for cosmologists! 
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